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Abstract. We analyze the effective potential for nanoparticles trapped at a fluid 
interface within a simple model which incorporates surface and line tensions as well as 
a thermal average over interface fluctuations (capillary waves) . For a single colloid, a 
reduced steepness of the potential well hindering movements out of the interface plane 
compared to rigid interface models is observed, and an instability of the capillary 
wave partition sum in case of negative line tensions is pointed out. For two colloids, 
averaging over the capillary waves leads to an effective Casimir-type interaction which 
is long-ranged, power-like in the inverse distance but whose power sensitively depends 
on possible restrictions of the colloid degress of freedom. A nonzero line tension leads 
to changes in the magnitude but not in the functional form of the effective potential 
asymptotics. 
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1. Introduction 



The effective forces between rigid objects immersed in a fluctuating medium have 
attracted a steadily growing interest because their understanding allows one to design 
and tune them by choosing suitable media and boundary conditions and by varying the 
thermodynamic state of the medium. Possible applications range from micromechanical 
systems to colloidal suspensions and embedded biological macromolecules. Accordingly, 
these fluctuations may be the zero-temperature, long-ranged quantum fluctuations of 
the electromagnetic fields giving rise to the original Casimir effect [1] between flat or 
corrugated immersed metallic bodies [21 [3]. Other examples for fluctuation induced 
long-ranged effective forces between immersed objects involve media such as bulk fluids 
near their critical point [HE], membranes [6] or interfaces [7]. 
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A quasi two-dimensional realization of a fluctuating system with long-ranged 
correlations is given by a liquid-liquid or liquid-vapor interface along with its capillary 
wave excitations, for an experimental study in real space see Ref. |H]. Consequently 
rigid objects such as colloids which are trapped at the interface are subject to an 
effective force generated by the fluctuating capillary waves. Theoretical investigations 
of this phenomenon have dealt with the effective interaction between pointlike |^ or 
rodlike [9] objects. If the particles are held fixed, and the interface is pinned at 
their surfaces, the capillary wave mediated interaction between the particles is a direct 
analogon of the original Casimir effect for a two-dimensional Gaussian scalar field with 
Dirichlet boundary conditions at the one-dimensional boundaries, given by the three- 
phase contact lines (for a general treatment on the Casimir effect in scalar fields, see 
Ref. [10]). However, the interface and the colloids are embedded in three-dimensional 
space, and therefore the colloids and the position of the three-phase contact lines may 
fluctuate. It has been shown in Refs. [HI [121 [13] that this situation corresponds 
theoretically to a Casimir problem with fluctuating boundary conditions. Depending 
on the type of admitted boundary fluctuations, the asymptotics of the resulting Casimir 
interaction varies considerably. This is an effect absent in the case of objects immersed in 
three-dimensional bulk systems. Moreover, the fluctuations of the colloids in direction 
perpendicular to the interface plane influence their stability, i.e. the activation energy 
to desorb them from the interface. This aspect has received little attention so far, as 
often rigid-interface models are employed to discuss the energetics of trapped colloids 
at interfaces (see, e.g., Refs. [T^ fT5]). 

In this work, both questions, the stability of single colloids and their Casimir-like 
interactions in the presence of fluctuating capillary waves, are analyzed with particular 
attention to the effects of a nonzero line tension. The paper is structured as follows: 
In Sec. |2] we discuss a free energy model for the entrapment and stability of a single 
colloid at a fluid interface. In order to be self-contained, the predictions of the rigid 
interface models incorporating surface and line tensions are summarized in Subsecs. 12.11 
and 12.21 whereas the effect of capillary waves is treated in Subsec. 12.31 In Sec. |3] we 
discuss the fluctuation-induced interaction between two trapped colloids, making use 
of the fluctuation Hamiltonians derived for the single colloid problem. The importance 
of the type of admitted boundary fluctuations on the effective interaction is highlighted 
with the exemplary discussion of two cases: (i) both colloids fluctuate freely and (ii) 
the colloids are held fixed. For both cases, line tensions modify the magnitude of the 
effective interaction but are shown not to lead to qualitatively new behavior. 

2. Free energy model for spherical colloids trapped at interfaces 

In this section, we consider the free energy of a single, spherical colloid with radius 
R trapped at an interface between two phases I and II which is assumed to be flat 
in equilibrium. The height of the center of the colloid above the interface is denoted 
by 2, and we look for the free energy F of the colloid with z fixed, i.e. F = F{z) 
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Figure 1. Side view of a single colloid with radius R trapped at an interface (z = 0). 
The full circle shows the colloid in equilibrium, with ist center located at a height Zoq 
above the interface. The radius of the circular three-phase contact line is given by 
^o,cq- The dashed circle shows the colloid out of equilibrium at a center position z. 
The surface tensions phase I/phase II, colloid/phase I, and colloid/phase II are denoted 
by 7, 7i, and 711, respectively. The physical contact angle 9 differs from Young's angle 
cos^o ~ (71 ~ 7ii)/7 foi' a non-vanishing line tension. 



describes a constrained free energy and the equilibrium free energy for the colloid follows 
upon minimization with respect to z. This constrained free energy can be viewed 
as an effective potential, in which the colloid fluctuates. In Subsecs. 12.11 and 12.21 we 
introduce the surface and line tension contributions to F{z) under the assumption that 
the interface stays rigidly fiat. In Subsec. 12.31 the effect of interface fluctuations on 
the constrained free energy F{z) is calculated in a perturbative manner, employing the 
capillary wave concept. 

2.1. Surface tension 

Let 7 denote the surface tension of the fluid interface and 7i[n] denote the surface tension 
of the colloid surface (with area ^ipij) exposed to phase I [II] (see Fig. [1]). The free 
energy F{z)., measured with respect to the configuration where the colloid is completely 
immersed in phase I, is given by [11] 



Here, A is the cross section area of the colloid with the interface. Introducing the reduced 
quantities z = z/R, F{z) = F{z) / [tt'^R^], and Young's angle Oq via cos 6^0 = (71 — 7ii)/7, 
the free energy becomes 



According to this model, the free energy (or single colloid effective potential) is harmonic 
with spring constant k = F"{zeq) = 2-717. The equilibrium position is given by 
Zeq = — COS Oq and the depth of the harmonic well defines the activation energy 
Fa. = vr7i?^(l — I cos 6*0!)^). Thus we see that for interfaces of simple fluids (7 ~ 10~^ 
N/m) the activation energy becomes comparable to the thermal energy = k-^T only 




z = 






(2) 
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for nanoscopic colloids with radius i? ~ 1 nm. At such small scales the application of the 
simple, thermodynamic concept of surface tension seems doubtful. Nonetheless, recent 
investigation using computer simulations [lEl |T7] have indicated that the behavior of 
truly nanoscale particles can be described within a wide range of conditions by such a 
phenomenological thermodynamic model upon introduction of particle size dependent 
surface tensions and three-phase line tensions. 

2.2. Line tension 

The line tension was introduced by Gibbs [HI [19] to define the excess free energy 
associated to the line where three phases meet. The accurate experimental measurement 
of the line tension has been a considerable challenge. In fact, the uncertainty in the order 
of magnitude of the line tension has generated a considerable number of studies. The 
interested reader is referred to a recent review devoted to the current status of the 
three-phase line tension [20]. On theoretical grounds the line tension can be either 
positive or negative and it is expected to be a small force [19], 10"^^ N. For simple fluids 
away from criticality, dimensional analysis gives r ~ k^T/a where o" ~ 0.1 nm is a 
typical atomic length scale. Line tensions inferred from experiments span several orders 
of magnitude, 10~^^ — 10~^ N [20} |2T]. reflecting the variety of experimental techniques 
(which are all indirect and where systematic errors can hardly be estimated) and the 
variety of materials used. More recent experimental efforts point towards line tensions 
with magnitudes at the lower end of the range given above (see, e.g., the discussion in 
Ref. [22] and references therein). It is not clear from the outset that the concept of line 
tension is unique, i.e. independent from notional shifts of the interfaces (which also shift 
the location of the contact line) that should be permitted within the molecularly diffuse 
interface region. This question is treated in detail in Ref. [22] which lines out a protocol 
how to define and relate line tensions obtained in different experimental setups. 

For our purposes, we may assume that a specific definition for the interface locations 
has been adopted with regard to which surface tensions are defined (see Fig. [1]). If we 
associate a line tension r with the energy of the three-phase contact line (which is a 
circle with radius tq = y/R^—z^), the free energy of Eq. ([2]) must be amended as follows: 




(4) 



The spring constant k = F"{zcq) of the potential well near the minimum is also modified 
and becomes 
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Figure 2. Plots of the reduced effective potential F{z) for 4 different values of the 
reduced line tension f . With increasing f , the configuration crosses from the stable 
via the metastable to the unstable regime. Young's angle is set to 6*0 = 120° which 
corresponds to the contact angle of /im-sized polystyrene colloids at an oil/ water 
interface 



met ast ability for the trapped colloid [TS]. Large positive line tensions lead to the 
desorption of the colloid from the interface (see Fig. [2] for the variation of the effective 
potential F{z) with f for 6o = 120°). Thus, within this simple model, a first order 
transition between colloids desorbed in one of the bulk phases and adsorption to the 
interface is possible [T5] . 

2.3. Capillary waves 

In thermal equilibrium, the interface is not sharp but acquires a finite thickness through 
density fluctuations. In a coarse-grained picture, these density fluctuations correspond 
to fluctuations of the mean interface position u(x) around the equilibrium position u = 0. 
Here, x = (x, y) is a vector in the equilibrium interface plane z = 0. For wavelengths of 
these interface position fluctuations which are larger than the correlation length in the 
bulk phases, the free energy of an interface configuration m(x) is given by the surface 
energy of this configuration: 



Here, S^cn is the interface or meniscus area projected onto the plane 2; = of the 
reference interface, i.e. it is the whole plane z = with colloids which are possibly 
trapped at the interface "cut out". The free energy in Eq. contains an additional 
term which accounts for the costs in gravitational energy associated with the meniscus 
fluctuations. It involves the capillary length Ac given by Ac = [7/(|pii — Pilfl')]^^^, where g 
is the gravitational constant and pi the mass density in phase i. Usually, in simple fluids 
Ac is in the range of millimeters. Surface excitations whose free energy are described by 
Eq. ([6]) are termed capillary waves, and their effect on static and dynamic properties of 




(6) 
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liquid interfaces is a subject of lively interest [21 [2^ [2^ [271 [251 [251 [3U] . 

The effective potential F{z) for the colloid whose center is fixed at height z above 
the equilibrium interface can be obtained through the partition function of the capillary 
waves: 

F{z)= -k^TlnZ, (7) 
Z = Zq^ [ Vu exp{-pn) ip-^ = k^T) . (8) 



Here, Zq is a suitable normalization factor. The Hamiltonian Ti. which enters the 
Boltzmann weight for a certain interface configuration u(x) is the difference in free 
energy between the configuration {u{x.), z} (describing the interface and colloid position) 
and the equilibrium configuration {u = 0, z^q} which we call the reference configuration. 
The equilibrium position of the colloid is determined through Eq. (j4|). Therefore 

H[m(x), z] = (F^en[{n(x), z}] - F^cn[{0, Zcq}]) + 7iA^i + 7iiA^ii + rAL (9) 

The difference in the interface areas colloid/phase I and colloid/phase II between 
the configuration {u, z} and the reference configuration is given by AAi and AAn, 
respectively, and the difference of the three-phase contact line length between these 
configurations is given by AL. In general, T-C[u{'x),z] is a complicated functional of 
u and z. In order to reduce it to a tractable expression which allows the analytical 
determination of the functional integral in Eq. ([8]), we perform a Taylor expansion in 
u and z [31]. To quadratic order, this implies that the Hamiltonian can be split into a 
two-dimensional "bulk" term Ticw and a one-dimensional "boundary" term Ti.-^,: 

7^[«(x), z] = Kw[«(x)]L,5_^^,^^ + 7^b[«(x), • (10) 

The two-dimensional "bulk" area S'men.ref is given by the interface in the reference 
configuration (i.e. the plane z = with the colloid cut out), and the boundary line 
c^5'men,ref is giveu by the three-phase contact line in the reference configuration. With 
these definitions Hew corresponds to the usual capillary wave Hamiltonian [32] : 

,2" 



7^cw[^^(x)] = I / d"^^ 

^ 'S'mcn.rcf 



:iii 



The boundary term Hh only depends on the difference h = z — Zeq and the vertical 
position of the contact line /. The latter is expressed by its Fourier transform: 

oo 

fi'f)=uidS„i)= ^-e^"^'^. (12) 

m=— oo 

In Eq. f[T2|) . the polar angle (p is defined on the reference contact line circle dSraen,rei- The 
Fourier coefficients Pm are referred to as contact line multipoles below, and since the 
contact line height f{ip) is real, Pm = Pl^n holds. With these definitions, the boundary 
term acquires the form 

7^b[M(x),/^] =7^b,i+^b,2, (13) 



2{Po-hy + Aj2\Pm\^ 



m=l 
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n 



b,2 



7r7 



T 



0,cq 



-2(Po-/^)' + 4^(m2-l) |P„ 



m=l 



(15) 



As before, ro,eq = \JR^ — -^cq is the radius of the circle enclosed by the reference contact 
line and f = t/{'-jR) is the reduced line tension. The detailed derivation of the two 



terms contributing to TYb can be found in [Appendix A[ Note that Tib i which describes 



the change in colloid surface energy upon shifting the contact line is strictly positive 
definite, whereas Hb,2 is not positive definite, regardless of the sign of r. 

The partition function Z (Eq. ([H])) can be written such that the integral over contact 
line fluctuations j Vf appears explicitly: 



Z = Z,^ j Vu exp (-/3KwN, z]) j Vf H 5[n(x) - /(x)] exp (-/3Hb[/, h]) . (16) 

X G iSni c n , r c f 

The integration measure for the contact line fluctuations is given by Vf = 
dPoY[m>od^^Pmdlio^P-m- However, in this form the 2d- "bulk" fluctuations u are not 
yet separated from the "boundary" fluctuations /. This can be achieved by splitting the 
field u of the local interface position into a mean-field and a fluctuation part, u = u^i+v. 
The mean-field part solves the Euler-Lagrange equation (—A + A~^) = with the 
boundary condition laSj^^^^^f = /. Consequently the fluctuation part vanishes at the 
contact line: v lawmen rot ~ 0- Then the partition function Z = ^fiuc^mf factorizes into a 
product of a fluctuation part independent of the boundary conditions and a mean field 
part which depends on the fluctuating contact line /: 



^fluc =Zq'Jvv Yl 5(t;(x))exp(-/?KwH) 

X^l95'|iinri.rcf 



(17) 



Z^i{h) 



Vfexp{-^f «/(x)(9„t<„,(x)) >.exp(-SHb|/,A]) , 



The first exponential in Z^f stems from applying Gauss' theorem to the energy 
associated with Umf- In this term 9„Umf denotes the normal derivative of the mean- 
field solution towards the interior of the circle i9S'men,ref; and di is the infinitesimal line 

segment on 5S'men,ref- 



Since Z^^c does not depend on the colloid position h = z 



"eq; 



it only contributes 



an additive constant to the effective potential. In Z^f, only the monopole fluctuations 
of the contact line Pq are coupled to h (see Eqs. ((Hj) and (fT5l) ). The mean-field energy 
term (the first exponential in Z^^, Eq. f|T8|) ) is diagonal in the multipole moments P^ 
(see [Appendix Bj) : 



2 fas. 



dif{x.) (a„Mmf(x)) ^"^^ -TTp-f 



p2 

Inl 



+ 2^m|P„ 



m>0 



with the reduced capillary length given by Ac = 1.12Ac/ro,cq- Therefore only the integral 
over the monopole fluctuation Pq yields a dependence on h in the partition function: 



Z = Zq[v, Pm{m^O)] / dPo exp <^ -TTp-f 



In Ac 



T 



0,cq 



(19) 
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exp < 



Pi 



In Ar 



0,eq 



-1/2 



X 



1 - 



r 



0,eq 



1 + In A J 1 



T 

0,eq 



(20) 



(21) 



The /i-independent contributions have been put into the new normahzation factor Zq. 
Thus the effective potential for the colloid moving around its equilibrium position is 
given by 

F{z) = -p-^ \nZ (22) 



1 - 



r 



const. + 7r7 (z — z, 



0,eq 



cq; 



In Xr 



T 

0,eq 



(23) 



The spring constant k = F"{zcq) of the effective potential can be compared to the one 
derived for the case of a rigid interface (Eq. ([5]). It is seen that we recover the latter upon 
neglecting the term involving the capillary length Ac- However, for a physical situation, 
the capillary length is not small and greatly diminishes the steepness of the potential 
well (for colloids with R = 10 nm at an air-water interface, — 90° and negligible 
line tensions, k is reduced by a factor of 14). In accordance with the Goldstone boson 
character of capillary waves, the spring constant vanishes for Ac ^ oo since the whole 
interface can move with no energy cost upon shifting the colloid. 

In accordance with the rigid interface result (Eq. ([5])), the effective potential 
becomes unstable for f > ^q. However, this perturbative calculation of the effective 
potential F{z) allows no conclusion on the height of the energy barrier, which exists 
in the metastable regime of positive line tensions (see Fig. [2]). One may speculate that 
the considerable reduction of the potential well steepness by capillary waves goes along 
with a reduction of the barrier height and would thus facilitate particle desorptionl|| 

For negative line tensions, the shape of the effective potential in Eq. ( l23l) is hardly 
affected since In Ac ^ 1. However, an instability shows up in the partition function 
which is contained in the factor Zq. This factor contains a contribution of the form 



Z'o oc V'fexp i -27r/?7 ^ \P„ 



m>0 



(m + 1) + (m^ - 1)- 



0,eq 



(24) 



where the measure for the contact line fluctuations does not contain the monopoles, 
Vf = J fl^^Q (iRePmC?IniP-m- Clearly, for arbitrarily small but negative r there exists 
a critical multipole order rUc above which the exponent becomes positive and thus the 
partition function becomes infinite. Taken at face value, for negative line tensions the 
interface would become unstable by forming ripples with small wavelengths near the 



I This effect might be expected quite in similarity to the effective reduction of the barrier height in a 
double- well potential for a single quantum-mechanical particle [33] . 
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colloid. In a physical situation, there is however a lower cutoff in the wavelength of 
these ripples set by the colloid surface roughness which presumably adds a positive 
energy penalty to higher multipole fluctuations of the contact line. 

As is well-known, the capillary length Ac serves effectively as an infrared cutoff 
for the capillary wave spectrum. Although the capillary wave model used in this work 
was derived assuming gravitational damping of the capillary waves (see Eq. ([6])), the 
dependence of all resulting expressions on Ac is not specific to that form of the capillary 
wave Hamiltonian [31]. Such an infrared cutoff results equally well from a finite system 
size (e.g. a two-phase system with trapped colloids in a container with extension L 
or colloids trapped on a droplet with radius -Rd)- Thus for such systems, the effective 
potential (Eq. (123|) ) is obtained by just replacing Ac by L or R^. Since for smaller system 
sizes the influence of the line tension on F{z) is more pronounced, it is conceivable to 
obtain a value for r from the fluctuations in the colloid position which sample F{z). As 
can be seen from Eq. (1231) . no knowledge of the surface tension 71 and 7n (or Young's 
angle 6*0) is required, but the modified contact angle | sin6'| = fo,eq enters. Certainly, an 
experimental realization appears to be difficult because of the difficulty in determining 6, 
and one could resort to simulations in a first step which determine the colloid fluctuations 
with varying system size and in which 6 can be determined straightforwardly. 

3. Fluctuation induced forces between two colloids 

The previous considerations can be extended to the case of two colloids which are 
trapped at the interface at distance d. Clearly, if both colloids are at their equilibrium 
position (defined by Eq. (jl])) and capillary waves are neglected, the interface is flat and 
therefore no interface-mediated interactions are present. If, by some external force, the 
colloids are moved away from equilibrium, the interface will adapt to a long-ranged 
deformation and induce a capillary interaction energy oc Inci between the colloids |34j . 
the well-known cheerios effect [35] . The occurence of such a long-ranged interaction is 
tied to the occurence of a net force on the system "colloids + interface" [36l [3ll 137] . 

Here, we are interested in the occurence of an interface-mediated interaction 
potential V{d) in the force-free situation which are brought about by the fluctuating 
capillary waves. To that end, we can apply the partition function analysis developed in 
the previous section, extended to the case of two spherical "obstacles" trapped within 
the interface. We will focus on two scenarios: 

(Al) No external force acts on the colloids, therefore the colloids are free to fluctuate in 
the direction perpendicular to the interface. 

(A2) The colloids are fixed at their equilibrium position by external means. On average, 
there is no external force acting vertically on the colloids, although at a given 
instant of time some force is needed to counteract the Brownian fluctuations of 
the colloids. In this case, the effective potential V{d) between the colloids is 
related to the pair correlation function g{zeq, Zeq,d) between the colloids through 
l3V{d) = -\ng{zeq,Zeq,d). 




Figure 3. (a) Side view of two colloids at distance d trapped at an interface, with 
hi {i = 1,2) denoting the relative center position of the fluctuating colloid i (dashed 
circles) with respect to the colloid i in equilibrium (full circles), (b) Top view on the 
equilibrium interface, spanning the x — y-plane. The colloids in equilibrium occupy 
the circular areas S^, and their boundaries dS^ arc the equilibrium (or reference) 
three-phase contact lines. 



In both scenarios, the interface and in particular the three-phase contact hne are free 
to fluctuate, subject to the energy penalty of the capillary wave and the boundary 
Hamiltonian lined out in the previous section. In previous work [121 IIS] (neglecting line 
tensions) we have established that the long-range behavior of the effective potential 
V{d) depends sensitively on the types of contact line fluctuations. For the case of a 
pinned contact line on the colloid and the colloids fixed at their equihbrium position, 
V{d) oc Inlnrf whereas for the case (Al), V{d) oc d~^. In the present work, we will show 
(i) that for case (A2), fixed colloids but unpinned contact line, the effective potential 
is still long-ranged, V{d) oc ln(l + Ind) and (ii) that line tensions do not change the 
leading power in the long-range behavior of V{d) but the instability for negative line 
tensions in the partition sum for a single colloid also occurs in the leading interaction 
term for case (Al). 

The effective potential V{d) is obtained via the partition function of the fluctuating 
capillary waves Z via 

V{d) = -k^T\nZ{d) . (25) 

The partition function Z oc J Vu exp(— /37i) contains a Hamiltonian which as before 
contains a 2d- "bulk" term and a sum of boundary terms for each of the two colloids: 

2 

nH^),z] = Kw[n(x)]|,,5_^^^^ + ^HUn(x),2;]|^^^^o ■ (26) 

i=l 

The functional form of Ticw and the 7i\^ is given by Eqs. f|TT]) and f|T3|) - f|T5|) . respectively, 
with due generalization of contact line multipoles for each colloid z, Pim, and of 
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colloid height differences h ^ hi = Zi — z^q- In the equilibrium (reference) configuration, 
colloid i intersects the interface plane in the circular area S^, thus the three-phase 
contact lines in the reference configuration are given by dSf {i = 1,2). Thus the 2d- 
"bulk" area over which the capillary waves fluctuate is given by S^cn,rcf = \ ULi 
(see also Fig. [3]for the geometric deflnitions). Via the integration domain Smon.rof of Ticw, 
the total Hamiltonian and hence the partition function Z[d) of the system depends on 
the distance d of the colloid centers. 

As before, the fluctuations over the contact lines fi = exp(im</)i)Pjm are 
incorporated into the partition function via (5-function constraints: 

Z{d) = Zq-i j Vu exp {-mc^u, d]) X (27) 

n/^/^ n 5[u{^^)-f^{^i)]^^v[-m^[f^M) ■ 

The normalisation factor is chosen such that Z{d oo) 1. The difference between 
cases (Al) and (A2) deflned above shows up in the deflnition of the measure for the 
contact line fluctuations: 

{dhidPioYl^^f^dRePimdlmPi_m (Al) 
(28) 
dPio Um>o dRePimdlmPi.m (A2) 

As seen above, in the unconstrained case (Al) an additional integral over the colloid 
height variables is performed. Both cases can be discussed conveniently by splitting the 
fleld u of the local interface position into a mean-fleld and a fluctuation part, u = u^^+v. 
The mean-fleld part solves the Euler-Lagrange equation (—A + A~^) = with the 
boundary condition u^i \dso = fi- Consequently the fluctuation part vanishes at the 
contact line: v Iq^o = 0. Then the partition function Z = Z^^^-Zjai factorises into a 
product of a fluctuation part independent of the boundary conditions and a mean fleld 
part which depends on the fluctuating boundary conditions fi of the meniscus on the 
colloid surfaces: 



^fluc = ^o"' /"pi; JJ JJ 5(i;(x,))exp(-/3Hcw[f,rf]) , 
i=i x,ea5f 

^mf = n / ^/^expj-^^ / ^rf£,/,(x,)(9„u^f(x,;d))|exp(-/57i;j,[/„/i,]) .(29) 

The flrst exponential in Z^i stems from applying Gauss' theorem to the energy 
associated with Umf- In this term dnUm{ denotes the normal derivative of the mean- 
fleld solution towards the interior of the circle dSi, and dii is the inflnitesimal line 
segment on dSi . 

In this form, the partition function is amenable to analytical expansions for small 
and large distances d between the colloids. The multiplicative separation of Z allows 
to deflne additive contributions to the effective potential: V = Vauc + Vmf with 



Stability and interactions of nanocolloids at fluid interfaces: effects of capillary waves and line tensions 12 




Figure 4. Pictorial representation of the separability of the total effective potential 
V caused by the fluctuating capillary waves and the possibly fluctuating colloid into 
a fluctuation (Vhuc) and a mean-field part (ymen)- The fluctuation part (second 
picture) is obtained by summing over all admissible capillary waves which are pinned 
at the colloid surface {u = 0) and with the colloids themselves fixed. The mean-field 
part (third picture) is obtained by summing over the colloid height fluctuations (long 
arrows) and the contact line fluctuations (short arrows) equipped with a Boltzmann 
factor incorporating the energy of the mean-field meniscus solution whose boundary 
condition is set by the momentary position of the contact line. 

/3 Vfluc[mf] = In 2^fiuc[mf] • The techniques to evaluate the fluctuation and mean-field part to 
the effective potential have been presented in detail in Ref. [13], and we give a summary 
of the main results which are necessary to discuss the influence of the line tension terms 
on V{d). 

3.1. Fluctuation part 

The fluctuation part contributes equally for both cases (Al) and (A2) introduced above. 
The 5-functions in the fluctuation part of the partition function can be removed by 
using their integral representation via auxihary fields ipi{^i) defined on the interface 
boundaries dS^ [38]. This enables us to integrate out the field u leading to 




We note in passing, that the fluctuation part in the form of Eq. (l30l) resembles 2d 
screened electrostatics: it is the partition function of a system of fluctuating charge 
densities ipi residing on the contact circles. For large d/r^ it can be calculated by 
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utilizing the multipole expansion 

oo 

i^ii^i) = ^ Vimexp(im(y9i) , (31) 

m=— oo 

introducing the auxihary multipole moments ipim of order m pertaining to colloid i. 
Using these it can be shown that in the limit Xc ^ d ^ R the fluctuation part of the 
effective potential /3Vauc = — In-Zfiuc has the form 

oo 

ar, 



pVfinc{d) = const. + ao(c^) + ^ (32) 

n=l 

ao{d) = ^Inln-^ , (33) 

^ ^ 0,eq 

where the a„ for n > are numerical coefficients. Through the multipole analysis it is 
found that interaction terms oc ipim'>p2m' contribute terms to l^uc which are proportional 
to c^-2(m+m')_ xhe combination oc ^/'ioV^2o (fluctuating auxiliary monopoles) gives rise to 
the leading term oq in Vauc- 

3.2. Mean- field part 

The calculation of Z^f (Eq. (l29l) ) requires to determine the solution of the differential 
equation 

(-A + \f) M„,f = (34) 
with the boundary conditions at the fluctuating contact line and at infinity, respectively: 

Mmf(Xi)|^^gg50 = fi{^i) (35) 
^imf(x)||^|^^ ^0. (36) 

We write the solution as a superposition = ui + U2 where Ui = 

'Ylim-^rriiji/ X^Ai^e^'^'^^ is the general mean-field solution in \ (see Fig. [3] for 
the geometric definitions). The solution has to match to the boundary conditions at 
both circles dS^ and 882. This can be achieved by a projection of U2 onto the complete 
set of functions on dS^, {e""'^^}, and vice versa. Equating this expansion with the 
contact line multipole expansion fi{(pi) = 'YliiPi exp(im(y9j) leads to a system of linear 
equations for the expansion coefficients {Ajm}. This system can be solved analytically 
within a systematic 1/d expansion or numerically, observing rapid convergence. Owing 
to the linearity of Eq. ( IMl) . the mean field part of the partition function can be 
written in a Gaussian form: 

2mf = j Vfi exp (-/3HKf, d]) exp |-7r/37 - t^- j Y.^P^^ ' ^*)'| ' (^7) 

where 7i[umf] is a symmetric quadratic form in the vector of the contact line multipole 
moments = (. . . , P^o, Pn, ■ ■ ■)■ 
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Using this form, it can be shown that the mean-field part of the effective potential 
has a similar expansion to the one of the fiuctuation part Vquc (Eq- fl32l) ) in the limit 
Ac > > i?: 

oo , 

pV^,{d) = const. + ho{d) + (39) 

n=l 

Also similar to the analysis of Vfjuc? interaction terms oc PimP2m' contribute terms to 
Vmf which are proportional to c^-sCm+m') rpj^^ fiuctuating contact line monopoles and 
the possibly (in case (Al)) fiuctuating colloid heights hi gives rise to the leading term 
60 in The form of bo{d) and the values of the numerical coefficients fe„ depend on 
the cases (Al) and (A2) introduced above. 

(Al) Here, for the freely fiuctuating colloid, the integration measure was given by 
Vfi = dhidii. Upon change of variables hi ^ hi — PiQ, it is seen from Eq. (1371) 
that the ci-dependent part of -Emf is given by det E. For the case of vanishing line 
tension, the properties of det E were discussed in detail in Ref . [13] . In particular 
it turns out that the four leading terms in the expansion of Knf and Vhuc cancel 
each other {bm = —am for m = 0, 1,2,3) and the total effective potential is given 
to leading order by 

PV{d) « . (40) 

This corresponds to a quadrupole-quadrupole interaction according to the power 
counting in terms of the multipoles of either the fluctuating auxiliary charge 
densities ipim (for Vauc) or the fluctuating contact fine Pim (for l^mf)- Hence, a 
nonvanishing line tension does not change the leading power in V{d) since the line 
tension contributions (via 7^6,2? see Eq. (1151) ) to detE are nonzero only for contact 
line multipoles higher than dipolesJ§| Therefore, including the line tension we find 
for the total effective potential 



3f 




3f 


0,cq 




' 0,cq 












)■ 




0,cq , 





pV{d) = const. '^^^^ 4 (41) 



The dependence of the effective potential on r is depicted in Fig. [51 For positive line 
tensions, l3V{d) is always asymptotically attractive, save for the value f/r^^^ = 1/3 
where the coefficient of the leading d~^ term vanishes, and the effective potential 
becomes even shorter ranged. For negative line tensions, we encounter a divergence 
of this coefficient for f/rgg^ — > —1. This divergence is related to the instability 
in the one-colloid partition function at the interface already discussed following 
Eq. (El. 

§ The boundary Hamiltonian contains also a line tension contribution for the monopole terms 
oc (Pjo — (see Eq. ([T5|) ). Upon integration over the colloid height /i^, the line tension dependence 
here is absorbed in a d-indepcndcnt multiplicative factor in Z,„f , see Eq. ([37]). 
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Figure 5. For case (Al), the freely fluctuating colloid, the dependence of the cocfScient 
of the leading term (04 + 64) in the effective potential V{d) on the line tension is 
shown. The dependence on r enters through the reduced variable x = t/tq^^ where 
f = t/{'jR) and fo,oq = ?'o,cq/i?- 



(A2) This case implies fixing the colloids at their equilibrium positions Zcq. Thus in the 
integration measure for the contact line fluctuations the integration over the colloid 
height h = z — Zcq is absent, Vfi = dti. According to Eq. ( 1371) . the quadratic form 
f^E f in the Gaussian integral is slightly changed by the second exponential on the 
right hand side of Eq. fl37|) . In particular, this leads to a changed leading coefficient 
ho in the mean-field free energy 



ho{d) 



- In In — 

2 ro 



In 



eq 







1„ ' \ 










'0,cq_ 


^0,cq ) 



(42) 



Therefore, for case (A2) the leading term in the total effective potential contains a 
very long-ranged leading term of the form 



pV{d) 



const. H — In 
2 



0,eq 



In- 



d 



^0, 



(43) 



cq 



which for d ^ ^o,cq slowly approaches the asymptotic form found for Vquc 
(Eq. fl33l) ). Thus for fixed colloid position in the interface their pair correlation 
function g{zeq, Zeq,d) = — exp{—pV{d)) oc — ln((i/ro,cq) contains a long-ranged 
piece dominated by the fiuctuating "bulk" capillary waves only. This has also been 
found in a study treating the colloids as point particles [7] and in an analytical study 
of the pair correlation function in phase-separating 2d and 3d lattice models [39], 
and the asymptotics of the colloid pair correlation function is the same as exhibited 
by the fluid pair correlation function in the interface region |10]. However, treating 
the finite size of the colloids correctly leads to sizeable corrections in the asymptotics 
of the effective pair potential (Eq. ( l43l) ) and gives a nonzero V{d) in the physically 
relevant case (Al) (it is zero in the limit of point colloids). 
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4. Conclusion 

In this paper, we liave studied tlie influence of capillary waves on the stability and 
interactions of colloids (with radius R) trapped at a fluid interface with surface tension 
7, with particular attention to the effects of a line tension r. Quite often, the stability 
of colloids at a fluid interface with respect to vertical displacements h from their 
equilibrium position is discussed using a rigid interface model. This gives for negligible 
line tensions and partially wetting colloids a steep potential well with spring constant 
k = 27i'-f. A finite line tension changes the spring constant by a term (x r/R and may 
induce metastability for the trapped colloids for certain positive values of r. Within a 
perturbative model we have found that the potential well is considerably broadened by 
capillary waves (qualitatively, k — > k/\n{\c/R) where Ac is the capillary length in the 
interface system). This suggests also a reduction of the metastability barriers in case of 
positive line tensions, although calculations beyond our perturbative model (quadratic 
in the fluctuations) are needed for conclusive results. 

Capillary waves also induce effective interactions between two colloids which are 
of Casimir type. For freely fluctuating colloids a power-law dependence of the effective 
potential V{d) in the intercolloidal distance d is obtained, V oc d~^. A flnite line tension 
does not change this power-law dependence, save for a speciflc positive value of r where 
the corresponding coefficient vanishes and V{d) becomes even shorter-ranged, decaying 
at least oc d~^^. Negative line tensions increase the amplitude of V{d). For colloids 
fixed in the interface, the effective potential is equivalent to the potential of mean force 
between them and acquires a long-ranged component V{d) oc ln(l + Alnd) where A is 
a line-tension dependent coefficient (see Eq. fH3l) ). For d ^ R, our results contain as 
a special case the long-ranged potential of mean force already discussed for pointlike 
colloids within an fluctuating interface. 

In previous work [121 [13] we have discussed the strong attractive component in the 
fluctuation force which occurs for small separations between the colloid. This strong 
attraction is independent of the surface properties and also of the line tension and can 
be understood from the capillary wave partition sum with stricht Dirichlet boundary 
conditions on the colloid surface (see Subsec. 13.11) . Both short-ranged and long-ranged 
regimes of the effective fluctuation potential should be important for the aggregation 
of nanocolloids at interfaces and compete with other effective interactions such as of 
electrostatic origin [CT 112] . 

Acknowledgment: M. O. thanks the organizers of CODEF II for their invitation 
and the German Science Foundation for flnancial support through the Collaborative 
Research Centre SFB-TR6 "Colloids in External Fields, project section D6-NWG. 

Appendix A. Derivation of the boundary Hamiltonian 

In this appendix we derive the boundary term TYb which describes free energy changes 
upon shifting the contact line (cf. the result in Eqs. ( fT3l) - (fT5l) of Sec. 12. 3p . According 
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to Eqs. (I91)- fllll) the boundary term is given by 

= 7iAAi + 7iiAAn + 7AAproj + rAL (A.l) 

and contains contributions associated with the difference in the interface areas 
colloid/phase I and colloid/phase II between the configuration {u, z} and the reference 
configuration {u = 0, ^eq} (given by AAj and AAn, respectively) and the difference of 
the three-phase contact line length between these configurations (given by AL). The 
term AAproj describes the change in area (with respect to the reference configuration) 
of meniscus u{r) projected onto the plane z = 0. 

If the three phase contact line is slowly varying without overhangs, the following 
geometric relation holds between its projection onto the plane z = (parametrized in 
polar coordinates by ro{ip)) and the contact line uo = u{ro{ip),(p) itself: 

ro(^)= [R'-[u{ro{^),v)-z]'Y^' (A.2) 

= Keq + 2 2eq [«(ro(<^), ^) - h] - (u(ro(<^), (p) - hf] ^''^ . 



In Eq. ( ]A.2p . ro,eq = \/R^ — is the radius of the circular reference (or equilibrium) 



contact line, h = z — z^c^i^ the deviation of the colloid center from its equilibrium height 
and M(ro((/?), v^) is the actual height of the contact line parametrized in terms of the 
polar angle (p. 

Because of An = inR"^ — Ai, AAi = —AAu holds for fluctuations of the colloid 
surface area in contact with fluid I and II, respectively. Then, the associated changes of 
the free energy can be written as 

7iAAi + 7n AAn = 7 cos 60 d(p dr 



JO Jr,,., Vl - r2/i?2 

l'2n 

-R7cos6'o / dip [u{r(j{(p), ip) — h] (A. 3) 







7i?cos^o , r. ,n2 7^cos^o ^ \ 2( ^ 21 
d^ [f -h] / d^ [ro((^) - rg^^^fj 



22eq Jo 22eq JO 

where we have applied Eq. ( ]A.2[) . Following Ref. [3l], in the last line we have 
approximated the actual height of the contact line by the meniscus height at the reference 
contact circle SSmcn.ref , i-e. u{rQ{ip),ip) ^ u(ro^eq,V^) = /(v^)- Correction terms to this 
approximation are at least of third order in u and / [31] . 

The free energy contribution associated with the change in projected meniscus area 
can be written as 

7A^p,oj = 7 / dip drr = - dip [rl^^ - rl{ip)] . (A.4) 

^0 Jro{<f) ^ Jo 

Combining Eqs. (1A.3P and (lA.4p . applying again Eq. (lA.2p and using relation (jlj) for 
the equilibrium position of the colloid {RcosOq/ z^^ = T/(7ro^cq) — 1) we find 

7iAyli + 7„AAn + 7AAp,oj = ^ dip [f - h]^ - ^ dip [f - h] . (A.5) 

Jo '"O.eq Jo 
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The free energy contribution related to the length fluctuations of the contact line 
is written as 

r.27r 

),cq 



rAL = T dip 
'o 

f27r 

dip 



\Jro{ipY + [d^ro{ip)Y + [d^u{ro{ip>),ip)f - tq,, 



r 



2r, 



0,cq Jo 



r, 



0,cq 0,eq 

Comparing Eqs. (]A.5P - (lA.6p . we find that the linear terms cancel out because the 
equilibrium position and rg^eq are determined by Eq. (jl]). Inserting the decomposition 
of f{ip) from Eq. f|T2|) and performing the integrals over ip finally leads to the form 
Tit, — "^b.i +'^b,2 given in Eqs. f|T^ - ffT^ for the the total boundary Hamiltonian, where 
the two contributions read 



^b,l = — 



m>l 



TX 



b,2 



0,cq 



-2(Po-/^)' + 4^(m2-l) |P„ 



m=l 



(A.7) 



(A.8) 



and describe changes in colloid surface energy and in line energy, respectively, upon 
shifting the three phase contact line. 



Appendix B. Derivation of the mean— field energy term in Eq. ( llSft 

Let (r, ip) be polar coordinates in the equilibrium interface plane 2 = where r = 
is the center of the circle enclosed by the reference contact line. The solution to the 
mean-field equation (—A + A~^)Mmf = with the boundary condition Wmf ('^o.eq, V^) = 
fi.^) = Em,^mexp(imv9) is given by 



u,^i{r,ip)= ^ exp(imv9) 



m=— oo 



^0, 



X 



eq 



(B.i; 



where Km is the modified Bessel function of the second kind and order m. For 
nanocolloids, Ac ^ R, so that one can use the approximation 

-ln(x/C) (C~1.12) (m = 0) 
(2m-2)!!x~'" (m > 0) 

which is valid for x <^1. Thus we find (with /' = df/dr) 



Kmix) 



(B.2) 



2-K 



2-K 



p2 



f^V'(-Wmf(^0,cq,</')Mmf(ro,cq,V')) - I W^^ \ / 

^0,cq \ Ac/Tccqy 



2^m|P„ 



(B.3) 



m=l 



which immediately leads to Eq. (118 
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